Normal state properties of spin-orbit coupled Fermi gases in the upper branch of 

energy spectrum 
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We investigate normal state properties of spin-orbit coupled Fermi gases with repulsive s-wave 
interaction, in the absence of molecule formation, i.e., in the so-called "upper branch". Within 
the framework of random phase approximation, we derive analytical expressions for the quasi- 
particle lifetime t s , the effective mass ml, , and the Green's function renormalization factor Z s in the 
presence of Rashba spin-orbit coupling. In contrast to spin-orbit coupled electron gas with Coulomb 
interaction, we show that the modifications are dependent on the Rashba band index s, and occur 
in the first order of the spin-orbit coupling strength. We also calculate experimental observable such 
as spectral weight, density of state and specific heat, which exhibit significant differences from their 
counterparts without spin-orbit coupling. We expect our microscopic calculations of these Fermi 
liquid parameters would have the immediate applications to the spin-orbit coupled Fermi gases in 
the upper branch of the energy spectrum. 

PACS numbers: 03.75.Ss, 05.30.Fk, 67.85.Lm 



I. INTRODUCTION 

Motivated by the recent success on the evidence of 
Stoner ferromagnetism for the repulsive Fermi gas in the 
upper branch of the energy spectrum [l| , there has been 
increasing interest in the nature of uncondensed Fermi 
gas (free of molecules) within the repulsively interacting 
regime 0], which naturedly becomes the well-controlled 
platform for simulating Landau Fermi liquid. Much of 
the interest in ultracold atomic gases comes from their 
amazing tunability. A wide range of atomic physics and 
quantum optics technology provides unprecedented ma- 
nipulation of a variety of intriguing quantum phenomena. 
Based on the Berry phase effect [3| and its non-Abelian 
generalization [J], Spielman's group in NIST has suc- 
cessfully generated a synthetic external Abelian or non- 
Abelian gauge potential coupled to neutral atoms. Re- 
cent experiments have realized the atomic 40 K [5( or 6 Li 
@ gases with spin-orbit coupling (SOC). These achieve- 
ments will open a whole new avenue in cold atom physics 

The effect of SOC in fermionic systems has conse- 
quently become an important issue in recent years, and 
attracts a great deal of attentions in ultracold Fermi 
gases. Most of existing works are devoted to the effect of 
SOC on the superfiuid state with negative s-wave scat- 
tering length [LU - tl3| and the Bardeen-Cooper-Schricffcr 
(BCS) to Bose-Einstein condensation (BEC) crossover 
[l3-[l7|- Furthermore, the SOC give rise to a variety of 
topological phases, such as the quantum spin Hall state 
and the topological superfiuid [18-21]. In addition to the 
study of SOC effect on these symmetry-breaking or topo- 
logical phases, the normal state contains various poten- 
tial instabilities and deserves attention. The considera- 
tion of SOC systems in the framework of Fermi liquid 
theory is therefore desirable. As is well known, Lan- 
dau's Fermi-liquid theory provides a phenomenological 



approach to describe the properties of strongly interact- 
ing fermions. Fermi liquid parameters characterizing the 
renormalized many-body effective interactions have to be 
determined through experimental results. 

The purpose of this paper is to study key normal-state 
properties of two dimensional (2D) Fermi gases with SOC 
in the repulsive regime — their quasi-particle lifetime, ef- 
fective mass and Green's function renormalization factor. 
Following previous studies of Landau's Fermi-liquid the- 
ory including SOC [13, Hj|, we are attempted to build 
a microscopic foundation of phenomenological parame- 
ters. Therefore within the framework of random phase 
approximation (RPA), we derive analytical expressions 
for the quasi-particle lifetime t s , the effective mass m*, 
and the Green's function renormalization factor Z s for 
a 2D Fermi gases with repulsive s-wave interaction in 
the presence of Rashba SOC. To make contact with cur- 
rent experiments directly, we also calculate experimental 
observable such as spectral weight, density of state and 
specific heat, and discuss their corresponding experimen- 
tal signatures. We shall show that the modifications are 
dependent on the Rashba band index s denoting the two 
directions of the eigenspinors of the Rashba Hamiltonian. 

The paper is organized as follows. The model Hamilto- 
nian and the renormalizations due to s-wave interaction 
in the presence of SOC is discussed in Sec. II. In Sec. 
Ill, Starting from the RPA self-energy of the SOC Fermi 
liquid, we derived all the analytical formula of the Fermi 
liquid parameters in presence of SOC. The experimen- 
tal observable quantities such as the spectral function, 
density of states and specific heat are calculated. Sec. 
IV is devoted to discuss the experimental measurements 
of these fundamental parameters and their correspond- 
ing experimental signatures. The comparisons with the 
ordinary Fermi liquid are also presented. 
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II. PERTUBATIVE THEORY OF 2D FERMI 
GASES WITH s-WAVE INTERACTION IN THE 
PRESENCE OF RASHBA SOC 



resonance, the strength of the re puls ive interaction can 
be tuned within a wide range |27H32j]. 



Model Hamiltonian 



B. Renormalizations due to the s-wave interaction. 



We consider a 2D spin- 1/2 fermionic system with 
Rashba-type SOC and s-wave interaction, described by 
the model Hamiltonian 



rl — Hq + Hi- 



The non-interacting part "Ho reads as, 
L 2m 



p2 

n a = Y^ c th h a (z x p) • er - /i]< 



(1) 



(2) 



where c p = (c P) <|-, c p ^) T , fi = kp/2m is the chemical po- 
tential, a represents SOC strength and kp is the Fermi 
momentum in the absence of Rashba-type SOC. The re- 
duced Plank constant h is taken as 1 in this paper. The 
non-interacting Hamiltonian Hq can be diagonalized in 
the helicity bases 



= ±( 1 

y/2 \ ise^M 



,s = ±l, 



(3) 



where <^>(k) = arctan(fc J( /fc :E ) and s is the helicity of 
Fermi surfaces, which represents that the in-plane spin 
is parallel or antiparallel to the momentum. The dis- 
persion relations for two helical branches are £k,s = 
(k 2 + 2s/cfl|k| — hp) /2m, where fc# = ma corresponds to 
the recoil momentum in experiments 0, H[ . The Fermi 
surfaces are given by £k,s = 0, which yields two Fermi 
momenta k s = nkp — sk^ with k = \f \ + 7 2 . We have 
defined the dimensionless SOC strength 7 = fen/fop. Re- 
cently, the experimental realization of the SOC degen- 
erated Fermi gases have been reported [1-0] ■ By apply- 
ing a pair of laser beams to the ultracold 40 K or 6 Li 
atoms trapped in a anisotropic harmony trap, the equal 
weight combination of the Rashba-type and Dresselhaus- 
type SOC is generated. Their elegant experiments are 
performed in the weakly repulsive regime, which could 
provide the possibilities to study the SOC degenerated 
Fermi gases in the normal state. 
The interacting part reads as 



Hi = 2 5 



cPkd 2 p(i 2 q 1 1 

J2^y> c k+ q ,t c p-q-4 c P^ Ck ^' 



(4) 



where g = 27riVa s /3v27rm(z, which is controlled by the 
s-wave scattering length a s . Here N is the total atom 
number, C, z = y/l/ muj z is the confinement scale of the 
atomic cloud in the z direction perpendicular to the 2D 
plane, and m is the mass for the ultracold atoms. Notice 
that the universal properties of the low-energy interac- 
tion among ultra-cold atoms depend only on the scatter- 
ing length a s [23 - |26| . With the technique of Feshbach 



Let's consider the problem in the helicity bases |k, s) 
The non-interacting Green's function is 



G°(k, w ) 



w - £k,s + «sgn(w)0 H 



(5) 



The Dyson's equation expresses the relation between 
the non-interacting and interacting Green's functions in 
terms of the self-energy S s as 



G.(k,w) = 



1 



(6) 



All the many-body physics is contained in the self-energy 
S s . The poles of interacting Green's function G s (k, ui) 
give the quasi-particle excitations, of which lifetime t s = 
l/r s can be obtained from the imaginary part of self- 
energy as 



r s (k) = -2Im£ s (k,£ k , s ). 



(7) 



The real part of the self-energy gives a modification of 
dispersion relations. At low temperature, the properties 
of the low energy excitation in the vicinity of the Fermi 
surface is essential. Thus we can expand the real part of 
self-energy to the first order of lu and |k| — k s as 

ReS s (k, lu) = Re£ s (k s , 0) + ud u ReE s (k s , u) |„ 

+ (k-k s )-V k ReS s (k,0)| ks . (8) 

The interacting Green's function now becomes 



G s (k,w) 



(9) 



where £j* is the modified energy dispersion. The Green's 
function acquires a renormalized factor 



Z s = 



1 



(10) 



where A s = 9 w RcS s (k s , lu)\q. The Fermi velocity Vk 3 
9££ s /dk\k B can be calculated as follows 



dk 



\\k\=k. = Z s (Kk F /m + dkBeZs(k,0)\k.)- (11) 



For non-interacting Fermi gas, the Fermi velocity is u° = 
nkp/m. Therefore we introduce the effective mass via 
s /9A;||k|=fc s = nkp/m*. The effective mass in terms 
of self-energy is 



— = i- ( 1 + ^-3 fc Re£f (k s ,0) 
m Z s \ KkF 



(12) 
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FIG. 1. Feynman diagrams for the self-energy of the SOC 
Fermi liquid in the presence of s-wave interaction. The Feyn- 
man rules are defined under the helicity bases. The label s 
and r denote the helicity index. The self-energy is calculated 
within the framework of RPA 331. 
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FIG. 2. (Color online), (a): The contours of integration in 
the complex uj plane for the self-energy given by Eq. (119[) . For 
intermediate states with the energy < £k- q ,r < w, the pole 
of the propagator falls in the first quadrant, (b): Schematic of 
the deformation of the contour into the imaginary axis. The 
self-energy given by Eq. (I19|l is equal to the integration along 
the path CI with an additional contribution of the residue as 
shown by the integral path C2. 



The Eqs. Q, dTUJ> and are our starting points 
of microscopic calculations of normal states properties, 
which embody the main properties of a quasi-particle in 
the Landau theory of Fermi liquids. 



III. SOC FERMI LIQUID PARAMETERS WITH 
REPULSIVE s-WAVE INTERACTION 



A. RPA Self-energy 

To investigate the renormalization effects in two 
Rashba energy bands separately, it is convenient to work 
in the helicity bases. The interacting part of the Hamil- 
tonian in the helicity bases is rewritten as 



/ ' ss ;rr 
k,p,q 



(k, P, q)Vk +q)a /¥>p_ q ,r>¥>p,r¥>k, S . (13) 



where the interaction vertex V ss '- rr i (k, p, q) = 
fl , / ss '(^k,6'k+q)/ rr '(6'p,6'p-q)- Due to the presence 
of SOC, the spin is locked to momentum. The interac- 
tion vertex acquires an overlap factor f sg i ((9k, #k+q) and 
f rr ' (0p> 0p-q)' which is defined by 



(14) 



where 8^ and 6 p are the azimythal angles of k and p re- 
spectively. Within the framework of RPA, the interaction 
vertex is modified as (see Fig. [T]) 



^(k,P,q, W ) 



e(q,w) 



/ ss ' (#k, #k+q)/rr' (Op, Op-q), 

(15) 

where the dielectric function e(q, w) = 1 + gx(Q, w ) and 
%(q, w) is the bare density-density susceptibility of non- 
interacting SOC Fermi gas. In the long wavelength and 
low frequency limit, the susceptibility can be carried out 



as 



ftex(z/) = 

Imx(t/) = 



m 
m 



e(M-i) 



7T 



=0(1 - |y|), 



(16) 



where y — muJ / Kkp\<^\- It can be seen from Eq. (fT6"|) that 
the susceptibility satisfies the following relations 



x*(q,^) = x(-q.-w) = x(q,-w). 



(17) 



It's important to notice that the bare susceptibility is real 
for > vp^/l +7 2 |q|. For \uj\ < v F yJl + 7 2 |q|, the 
bare susceptibility in Eq. (|16l) contains an imaginary part 
which represents the absorptive behavior of the medium. 
This imaginary part is responsible for the finite lifetime of 
the quasi-particle in the medium. Along the imaginary 
axis, the analytical formula of x(q, uS) is much simpler 
(see Fig. [2]) as 



m . \y\ 

7T ^Jy 2 + 1 

The RPA self-energy (see Fig. [TJ is 



(18) 



S,(k,w) 



c 



d 2 qdqo 
(2tt) 3 



E 



e(q,?o) 



G°(k-q,u}-q ), 



where the J- sr is the overlap factor 

_ 1 + -srcos (6> k - 6> k _ q ) 

J~ sr — ^ 



(19) 



(20) 



The integral path C is shown in Fig. [2] (a). After defor- 
mation of contour path in Fig. [5] (b), the pole gives rise 
to a residue contribution 



Er e (k,^) 



E 



d 2 q g 



v W 



(21) 
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where the region of integration D r is k r < |k — q| < 
|k| - (r - s)k R and e = 1 + gxo (q, w - £k-q,r)- The line 
integral along the imaginary axis is 



j* r (k, U )=- rip 



d 2 qdq ^ gF s: 



1 



(2tt)3 ^ e ( qj i qo ) uj - iq Q - f k _ qr ' 
The total self-energy is given by 



line 



(23) 



The advantage of this decomposition is that the imagi- 
nary part is given by the contribution of residue, and the 
real part is mainly determined by the line integral. 



B. Quasi- particle lifetime 

The quasi-particle lifetime can be calculated by the 
imaginary part of self-energy, which comes from the con- 
tribution of residue in Eq. (|2ip. At zero temperature 
and for quasi-particle £ kjS > 0, the imaginary part of the 
self-energy reads 



r s (k) =2^e(6 c - q ,r)e(£ k)S - £ k _ q , r ) 



x IrnV^ r A (k, k - q, q; £ k , s - &_ q , r ). (24) 
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FIG. 3. (Color online). The inverse of the lifetime t s for 
40 K ultracold atoms as a function of the momentum k in the 
vicinity of the Fermi surface. The lifetime of quasi-particle is 
enhanced due to the presence of SOC. The parameters taken 
here are: the number of atoms is about 10 4 , ka = h/\ with 
A = 773nm, 7 = 0.5, trap frequency u z = 2n x 400Hz, and 
a a = 32ao, where ao is the Bohr Radius. 



wave repulsive interaction, where we can conclude that 
the quasi-particle is much stabler in the presence of SOC. 



The imaginary part of the RPA vertex in the helicity 
bases is 

ImV™? (k, k - q, q, k; w) = gT sr Im-}- , (25) 

where w = £ k , s — £ k -q.r and y — w/vp\q\. Since the main 
contribution of integral comes from the forward scatter- 
ing, i.e., y <C 1, the density-density susceptibility can be 
expanded about y = 0. One finds that the susceptibility 
in this neighborhood is 



x (y) = -(l + iy) + 0(y 2 ). 
The imaginary part of RPA vertex goes to 



(26) 



lmV s R r p r f (k, k - q, q; w) ~ - -y F sr . (27) 

Substituting Eq. (|2"7|) into Eq. ([2^1. the inverse lifetime 
T s (k) can be evaluated as 



r s (k) = - 



m 2 g 2 e F #(.5 1 o, 7 



, SUln 7 In-?, (28) 

ir(mg + ir) 2 \ 8 2 ' 4' 



where 5 = (k — k s )/kp. The result for ordinary Fermi 
liquid in the presence of s-wave repulsive interaction can 
be obtained by taking the limit 7 — > in Eq. (|2"5]l. The 
effect of SOC is shown in Fig. [3] via the comparison with 
the ordinary Fermi liquid with the same strength of s- 



C. Green's function renormalization factor 

Our starting point is the real part of the self-energy 
which contains two parts: one is the residue contribution 
given by Eq. ([2~T|). the other is the integral along the 
imaginary axis given by Eq. (|22[) . Thus we want to 
evaluate 

A s = A polc + A line 

= a ? sr 10 (k s ,OI?=o + a ? Ei inc (k s ,ok=o. (29) 

Given the rotation symmetry, we only need to consider 
k s = k s e x . The first term in Eq. (|2U)) is 



k — q,r) *F sr- 



(30) 



The second term in Eq. (|2"9")l can be integrated by parts, 
which gives two parts: the first one reads as 



1 line 



E/(0^-^)f^, (31) 



which cancels the residue contribution given by Eq. (|30D . 
Thus the final result of Eq. (|29|) can be expressed as 



A s = 



mgk 
(2" 



nfk 2 J d< ^ j d ^ J dx ^2 lm f sr ( x >y> 



(32) 
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FIG. 4. (Color online). Renormalization factor Z s as func- 
tions of scattering length a s with the same parameters in Fig. 
|3l The black thick line represents the renormalization factor 
of ordinary Fermi gases (7 = 0). The red (dashed) and blue 
(dotted) lines represent the analytical results for the SOC 
Fermi gas given by Eq. (|34p . The discrete points are evalu- 
ated numerically. 



where we have defined y = mw/qk s ,x — q/2k s and 

J sr 1 de 



iy - fj,(x, <j>) e 2 dy 



(33) 



1/2 + 



qir /|q|A: s , and the overlap fac- 
sr(l — 2x cos(j))/2l, with I = 



Here /z SjT .(a:, <j)) = 

tor F sr (x,(t)) = 

|k - q\/k s = ^/l - 4xcos0 + Ax 2 . 

From Eqs. (fTU)) , (|3"21) and ([5"3")l . we can obtain the renor- 
malization factor straightforwardly 



1 



1 



m z g z 



1 



8tt irng + ir) (k— S7) 2 



9 9 
m z g z 



8tt (mg + tt) 



(l + S7 + 0( 7 2 )). (34) 



The renormalization factor turns out to be dependent on 
the helicity s, and the leading correction due to SOC is 
0(7), which is different from the results of two dimen- 
sional electron gas (2DEG) with SOC in semiconductors 
[34j . We show the analytical results along with the nu- 
merical calculation for Z s as a function of the s-wave 
scattering length in Fig. 2] and the strength of SOC in 
Fig. [5] The Green's function renormalization factor for 
ordinary Fermi liquid is also shown for comparison in Fig. 
21 We can see from Fig. [S]that the renormalization fac- 
tor is reduced for the s — +1 branch while enhanced for 
the s = —1 branch with increasing strength of SOC. 



FIG. 5. (Color online). Renormalization factor Z s as func- 
tions of dimensionless SOC strength 7 with the same param- 
eters in Fig. [3] The red (dashed) and blue (dotted) lines 
represent the analytical results for the SOC Fermi gas given 
by Eq. (f34)l . The discrete points are evaluated numerically. 



isotropic due to the rotation symmetry. Without loss of 
generalities, we assume k s — k s e x in the following. We 
begin with 



r°° d 2 adw 
d k ReH s (k,0)\ lkl=k3 =ReJ 1 



(2tt)3 



yRPA 
v sr:sr 



(k, q, iw). (35) 



The interaction vertex is dependent on the external mo- 
mentum k because of the overlap of the helical eigen- 
states. It is instructive to consider some special cases. 
For weak SOC (7 <C 1), the integration in Eq. (|35"T) can 
be expanded to 

9 fe ReS s (k, 0) |, k|=fc , = + G(7 2 )- (36) 

In contrast to the SOC Fermi liquid with Coulomb inter- 
action [34:], our result is band dependent and has a first 
order correction 7. The effective mass reads as 



z-\i 



mg 
Attk 



,s 7 + 0(7 2 ))- 



(37) 



For strong SOC (7 ~ 1), we show the numerical results 
along with the analytical results of Eq. ([37)) in Fig. [6] 
The effect of SOC is shown in Fig. [71 from which we 
can see that the effective mass for the s = +1 branch 
is enhanced while the s = — 1 branch is reduced with 
increasing strength of SOC. 



D. Effective mass 



The effective mass can be evaluated by the real part 
of the static self-energy in the vicinity of the Fermi sur- 
face from Eq. (fT2")) . In contrast to the calculation of the 
renormalization factor Z s , the contribution of residue is 
irrelevant now. The correction of the effective mass is 



E. Spectral function, density of state and specific 
heat 



A close related quantity is the spectral function 
A(k, cj), which is the imaginary part of the single par- 



10 20 30 40 



0.25 0.5 0.75 1.0 
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FIG. 6. (Color online). Enhancements of effective mass as 
functions of s-wave scattering length a s with the same pa- 
rameters in Fig. [3] The red (dashed) and blue (dotted) line 
denote the case for s = +1 and s = — 1 respectively. The 
solid triangle and circle points are the corresponding numer- 
ical results. The black line represents the analytical results 
for ordinary Fermi liquid (7 = 0). The enhancements of the 
effective mass are dependent on the helical bands. 



FIG. 7. (Color online). Effective mass as functions of the 
strength of SOC 7 with the same parameters in Fig. [3] The 
red (dashed) and blue (dotted) line denote the case for s = +1 
and s = — 1 respectively. The solid triangle and circle points 
are the corresponding numerical results. The effective mass 
for s = +1 branch is enhanced while the s = — 1 branch is 
reduced with increasing strength of SOC. 



tide Green's function [371 \3 



A(k,w) 



-ImG ret (k,w), 



(38) 



where G ret (k, uj) is the retarded Green's function. It 
can be straightforwardly evaluated from the time-ordered 
Green's function G(k, ui) as 



ImG ret (k,w) 
ReG ret (k,w) 



ImG(k, uj)sign(uj) 
ReG(k,w). 



(39) 



With the great improvements in the spectroscopic tech- 
nique, it has been possible to directly measure the 
low-energy spectral weight function of a 2D system in 
the momentum-resolved radio frequency (rf ) experiments 
[4TJ - I42I ] . In Fig. |5J we show the spectral functions of the 
two Rashba bands separately. The spectral functions at 
the Fermi surfaces as shown in Fig. [5] (a) and (b) (vertical 
arrow) have the form 



A(k s 



(40) 



Fig. [8] (c) and (d) are density plots of the spectral func- 
tions, which could be compared with the results of the 
momentum-resolved rf spectroscopy in current experi- 
ments. 

Given the spectral functions, we can obtain the density 
of states (DOS) of the fermionic system with SOC via 

MM 



p(oj) = ^2p s (u;) = J2 J ■j^A.fru). (41) 



For the non-interacting case, the DOS now becomes 
'(J. 



p» = 



oj < - 



n y/2muj/k 2 F +K 2 



2m 



2m ' 
< UJ < 



' 2m ' 



(42) 



w > 



' 2m ' 



The DOS is modified in the presence of s-wave in- 
teraction. On the Fermi surfaces, for non-interacting 
SOC Fermi gas, the DOS in units of m/w is given by: 
p° +1 (E F ) = 0.28, pIxiEp) = 0.72, and p°(E F ) = 1 
for 7 = 0.5. For the interacting SOC Fermi gas, the 
DOS is evaluated numerically as follow: p + i(Ep) = 1.37, 
P-i(Ep) = 0.83 and p{E F ) = 2.20, where all the param- 
eters are the same with Fig. [3l 

The quasi-particles from two helical bands both con- 
tribute to the specific heat. At low temperature, the 
specific heat is proportional to the DOS on the Fermi 
surfaces and the temperature T. Thus the ratio between 
the specific heats at low temperature is 



P°(E F ) 



= 2.20, 



(43) 



where c° is the specific heat of SOC Fermi gases without 
interactions. 



IV. DISCUSSIONS AND SUMMARIES 

We have obtained various normal state quantities of 2D 
SOC Fermi gases in the presence of repulsive s-wave in- 
teraction. Ultimately, to make contact with experiments, 
two practical considerations warrant mention, (i): The 
repulsive s-wave interaction can be achieved on the up- 
per branch of a Feshbach resonance. One problem should 
be considered is that the upper branch of a Feshbach 
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FIG. 8. (Color online). Zero temperature spectral func- 
tion at different values of (k — k s )/kF are shown at (a) and 
(b). All the parameters are the same with Fig. [3] (a) 
is for s — +1 and (b) is for s = — 1 respectively. The 
eight peaks, from left to right, correspond to (k — k a )/kF = 
-0.01, -0.075, -0.05, -0.025, 0.025, 0.05, 0.075, 0.01. The 
vertical arrows in (a) and (b) denote S functions at uj = 
with weights 0.23 and 0.67 respectively, (c) and (d) are den- 
sity plots of the spectral function for the same parameters 
above, which correspond to s = +1 and s = — 1 respectively. 
The green line denotes the modified single particle dispersion. 



resonance is an excited branch, and will decay to the 
BEC molecule state due to inelastic three-body collisions 
[35l ]. However, with small scattering length, the decay 
rate is well suppressed P, Hi| and the system may be 
metastable for observation, (ii): Recently, the SOC de- 
generate Fermi gases have been realized in the ultra-cold 
atom systems [5|,|6j]. The effective SOC generated in their 
experimental schemes is an equal weight combination of 
Rashba-type and Dresselhaus-type SOC. In this paper, 
we have investigated the case of Rashba-type SOC. The 
Dresselhaus-type SOC is presented in Appendix, which 
is demonstrated to give the same normal state properties 
as the Rashba case. 

In current experiments, we consider the following typ- 
ical experimental parameters for quasi-2D systems. One 
can trap about 10 4 40 K atoms within a pancake-shaped 
harmonic potential with the trap frequency chosen as 
2tt x (10, 10, 400)Hz along the (x, y, z) direction. The 
system size can be estimated as (37.8, 37.8, 5.98)/im. 
The other related parameters are taken as: a s — 32ao, 
7 = 0.5, k R = 2n/\ = 8.128 x 10 6 m. The Fermi liquid 
parameters are listed in Table. U under this typical ex- 
perimental setup. For comparison, we give the results of 
ordinary Fermi liquid and the 2DEG in semiconductors 
together. We should notice that the single particle spec- 
tral function captures the valuable information for low 
energy excitations of the Fermi liquid and can be mea- 
sured by means of mometnum-resolved rf experiments 



[39M42j . The rf spectroscopy is a technique used to probe 
atomic correlation by exciting atoms from occupied hy- 
pcrfine states to another (usually empty) reference hy- 
pcrfine state. The single-particle spectral function is ob- 
tained in experiment through the mometnum-resolved rf- 
transfer strength. As a result, the inferences drawn from 
the spectral function in the vicinity of the Fermi surfaces 
can be used to determine the Fermi liquid parameters de- 
scribing the low energy behaviors of the normal state. Up 
to now, the rf experiment is the most promising method 
in ultracold atomic gases to measure these Fermi liquid 
parameters. 



TABLE I. Normal state properties for SOC Fermi liquid (7 = 
0.5), ordinary Fermi liquid (7 = 0), and 2DEG (7 = 0.051) in 
semiconductor. All other parameters used here are the same 
with Fig. El 



7 



0.5 



7 



2DEG( 7 = 0.051)* 



tt 



m*/m 



0.034kHz 0.037kHz 



Z+i = 0.23 
Z-i = 0.67 

m* +l /m 



4.88 
1.16 



0.55 



2.30 



55.36GHz 



0.97 



0.98 



f The results for 2DEG in InGaAs are taken from Ref. 
[34l | . Compared with the SOC Fermi liquid in the pres- 
ence of s-wave interaction, the results for the 2DEG 
with Coulomb interaction are independent on the en- 
ergy band, and the leading order correction relative to 
the SOC strength is 7 2 . 

ft The values of the inverse lifetime are evaluated at (k — 
k±i)/kF = 0.01. Compared with 2DEG in the typical 
semiconductor, the quasi-particle in ultracold atomic 
gases is much stabler. 



In summary, we studied the normal state properties 
of the SOC Fermi gas with repulsive s-wave interaction. 
The quasi-particle lifetime t s , the renormalization factor 
Z s , the effective mass m*/m are calculated, which em- 
body the main properties of a quasi-particle. To make 
contact with experiments directly, we calculated the spec- 
tral function A(k, lu), density of state p(w), and the spe- 
cific heat c v at low temperature. These quantities provide 
a good description of the low energy physics with SOC 
and s-wave interaction, which are measurable in current 
experiments. The analytical and numerical results show 
that the normal state properties are distinct for the two 
energy bands, and the leading correction relative to the 
ordinary Fermi liquid with s-wave interaction is on the 
order of 7, which are strikingly different from the SOC 
Fermi liquid with Coulomb interaction (34J. We expect 
our microscopic calculations of the Fermi liquid param- 
eters and related quantities would have the immediate 
applicability to the SOC Fermi gases in the upper branch 
of the energy spectrum. 
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Appendix : Relationships between the Rashba SOC 
and Dresselhaus SOC 



In this paper, the normal state properties of the Fermi 
liquid with Rashba-type SOC is studied. Now, we would 
like to demonstrate that another type of SOC, namely 
the Dresselhaus-type, gives exactly the same results for 
the normal state properties considered here. 

We start with the single particle Hamiltonian with 
Dresselhaus SOC 



Hd = 7^ + a{-p y a x - p x cr y ) - fi. 
The helicity bases with Dresselhaus SOC are 



\P,s) D 



~V2 



±1, 



ise 



(1) 



(2) 



where the subscript D represents the Dresselhaus-type 
SOC. The starting point of the microscopic calculation 
is the self-energy E s for each band. To illustrate the 
relationships of the two types of SOC, it is essential to 
derive the relationships of the Feynman rules between 
the two cases. For the Dresselhaus-type SOC, the non- 
interacting Green's function and the interaction vertex 



in the helicity bases reads 
G? fk,w 



w-£k, s +isgn(u;)0+' (3) 

v3. rr ,(K p, q) = gC Ok. ^+«)0 (V V*)- ( 4 ) 

The energy spectrum of the Dresselhaus-type SOC is the 
same with the Rashba-type. Thus the single particle re- 
tarded Green's function within the Dresselhaus represen- 
tation is the same with the result for Rashba-type SOC 
given in Eq. The overlap factor /^/(#k,#k+q) for 

the Dresselhaus-type SOC is given by 



f 



, q ) = i(l + ,s.sV^ (k) -^ k+q)1 ). (5) 



Compared with Eq. (fl"4"]l . we find the overlap factors are 
conjugated for the two cases 



ffi (#k, #k+q) — ffj (#k, #k-| 



(6) 



The bare susceptibility in Matsubara formalism is 
given by 



X 



D 



(q,iw n ) = k B T ^ G ? (k,iw m )Gf (k - q,iuj m -iuj n )T, 



D 

sr ' 



where the factor 



^ D 1 + sr cos ( 



(7) 



(8) 



is the same with the Rashba case in Eq. (|2"0)) . So the 
bare susceptibility x D f° r the Dresselhaus type SOC is 
equal to the case for the Rashba type SOC. The self- 
energy is given by Eq. (f2Tj) and (|22|) . The integrand 
function includes the following factors: the single parti- 
cle's Green's function, the bare susceptibility x°(qMi) 
and the overlap factor J 7 ® , which are all demonstrated to 
be the same for the two types of SOC. Therefore we con- 
clude that the normal state properties calculated here, 
such as the quasi-particle lifetime T s (k), the renormal- 
ization factor Z s and the effective mass m* /m are all the 
same exactly for the two types of SOC. 
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